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Recently [Europhys. Lett. 98, 37006 (2012)], based on heuristic arguments, it was conjectured 
that an intimate relation exists between the eigenfunction multifractal dimensions Dq of the eigen- 
states of critical random matrix ensembles D^i « qDq [q' + {q — q')Dq]~^ , 1 < g < 2. Here, we verify 
this relation by extensive numerical calculations on critical random matrix ensembles and extend its 
applicability to q < 1/2 and also to deterministic models producing multifractal eigenstates. We also 
demonstrate, for the scattering version of the power-law banded random matrix model at criticality, 
that the scaling exponents cr, of the inverse moments of Wigner delay times, (''w'') oc N~'^'' where 
A*' is the linear size of the system, are related to the level compressibility x as crq ~ 9(l^x)[l + 9X]~^ 
for a limited range of q; thus providing a way to probe level correlations by means of scattering 
experiments. 

PACS numbers: 05.45.Df, 71.30.-|-h, 73.23.-b 



I. INTRODUCTION 

It is well-known that the spatial fluctuations of the 
eigenstates in a disordered system at the Anderson- 
transition show multifractal characteristicsi*^ which has 
been demonstrated recently in a series of experiments)^ 
Therefore the modeling and analysis of multifractal states 
has become of central importance producing many inter- 
esting results. For this purpose random matrix models 
have been invoked and studied recently.—"— 

Since the exact, analytical prediction of the multi- 
fractal dimensions of the states for the experimentally 
relevant Anderson-transition in d = 3 or the integer 
quantum-Hall transition in d = 2 seems to be out of 
reach, it is desirable to search for heuristic relations in 
order to understand the complexity of the states at crit- 
icality. In the present paper wc propose such heuristic 
relations that are numerically verified using various en- 
sembles of random matrices. Moreover, we will demon- 
strate the robustness of these findings by extending to 
many different observables. These relations could play 
an important role in testing analytical predictions and 
could also be tested using the experiments mentioned 
above. 

The spatial fluctuations of the eigenstates can be de- 
scribed by a set of multifractal dimensions Dq defined by 
the scaling of the inverse mean eigenfunction participa- 
tion numbers as a function of system size N: 
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(1) 



where (• • •) is the average over some eigenvalue window 
and over random realizations of the matrix. For strongly 
localized eigenstates these quantities do not scale with 
system size, i.e. Dq for all q, while extended states 
always feel the entire system, i.e. Dq — )• d, for all q. Mul- 



tifractal states, on the other hand, should be described 
by the series of the _Dq, which are a nonlinear function of 
the parameter q. 

One of the most important generalized dimensions of- 
ten used in this context is the information dimension Di . 
It is defined through the scaling of the mean eigenfunc- 
tion entropy with the logarithm of the system size: 
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A further, well-known and widely used dimension is 
called the correlation dimension D2, which is extracted 
from the inverse participation number from Eq. ([T]) us- 
ing q = 2. The correlation dimension D2 shows up in 
the power-law scaling of the density-density correlation 
function in the energy domain as^ 
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as well as in the auto-correlation in space as^ 

d-D2 
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where and r; denote the position of sites k and /, 
respectively. In Eq. ((S]) the energy scale Eq is of the 
order of the bandwidthi^ 

Besides the wave functions the spectra also show 
anomalous properties. Spectral fluctuations can be char- 
acterized in many ways. A usual, often employed quan- 
tity is the level compressibility x, which is extracted from 
the asymptotic behavior of the spectral number variance 
as 
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where n{E) stands for the number of states in an inter- 
val of length E ii E ^ 1. The spectral fluctuations in 
a metallic system with extended states yield a vanishing 
compressibility, x ~^ 0, while in a strongly disordered 
insulating system the levels are uncorrelated, hence they 
are easily compressible, x = 1- However, for the multi- 
fractal states an intermediate statistics exists, < x < 1, 
furthermore the spectral and eigenstatc statistics are sup- 
posed to be coupled, which has been pointed out first in 
Ref. 

In a recent work^ Bogomolny and Giraud have shown 
that in a d-dimensional critical system the information 
dimension Di and the level compressibility x ^re related 
in a simple way as 



(6) 



furthermore the generalized dimensions Dq can be ex- 
pressed as 



r(fl-i/2). 
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X< 1 



These expressions have been shown to be valid for various 
critical random matrix ensembles in Ref. 

As for the critical, three-dimensional Anderson transi- 
tion and the two-dimensional quantum-Hall transition it 
has been shown earlier that there exists another relation 
between the level compressibility x the correlation 
dimension 



2x + D^/d = 1 . 



(8) 



This relation should obviously hold approximately only 
since the range of the generalized dimension and that of 
the level compressibility arc limited as < D2 / d < 1 but 
< X ^ li leaving the range of validity for the limit of 
weak-multifractality. 

Looking at Eqs. ([6|) and ^ it seems evident that there 
may exist further relations between generalized dimen- 
sions and the level compressibility, moreover, one may 
conjecture that the generalized dimensions of the multi- 
fractal spectrum must be also intimately linked together. 

In the present work we show a series of relations be- 
tween various generalized dimensions, and Dqi, and 
the level compressibility x allowing for a generalization 
that for particular cases yields Eq. © exactly and Eq. ^ 
in the appropriate limit. In order to prove that, numeri- 
cal simulations of various critical random matrix ensem- 
bles will be used. Moreover, by the use of our relations 
between generalized dimensions we state a clear link be- 
tween the spectral and scattering properties of disordered 
systems at the metal-insulator transition. Finally, by ex- 
ploring a deterministic model having a self-similar poten- 
tial that produces multifractal cigcnstates, we also show 
that our results arc not restricted to random matrix mod- 
els. 



II. MODEL AND HEURISTIC RELATIONS 

In Ref. Eqs. (gl) and © were shown to be correct 
numerically for the Power-Law Banded Random Matrix 
(PERM) model^iiiii^ at criticality. Below we will make 
use of this model to derive our main results. 

The PBRM model describes one-dimensional (Id) 
samples of length N with random long-range hoppings. 
This model is represented by iV x real symmetric 
(/? = 1) or complex hermitian (/3 = 2) matrices whose 
elements are statistically independent random variables 
drawn from a normal distribution with zero mean and a 
variance given by (|-ffmmP) = and 

1 1 



{\H„ 



(9) 



2 1 + [sin {n\m - n\/N) /{-Kh/N)]^^ 

where h and ^ are parameters. In Eq. (jQ]) the PBRM 
model is in its periodic version; i.e. the Id sample 
is in a ring geometry. Theoretical consider ation o^ ' ^ 
and detailed numerical investigations^iiiii^ have verified 
that the PBRM model undergoes a transition at /z = 1 
from localized states for /i > 1 to dclocalized states 
for /i < 1. This transition shows key features of the 
disorder driven Anderson metal-insulator transition;^ in- 
cluding multifractality of cigcnfunctions and non-trivial 
spectral statistics. Thus the PBRM model possesses a 
line of critical points h S (0,oo) in the case of /i = 1. 
In the following we will focus on the PBRM model at 
criticality, /_* = !. By tuning the parameter h, from 
h 1 to 6 ^ 1, the states cross over from the na- 
ture of strong- multifractality [Dq — > 0) which corre- 
sponds to localized-like or insulator-like states to weak- 
multifractality {Dq — >■ 1) showing rather extended, i.e. 
metallic-like states. Meanwhile at the true, Anderson 
transition in d = 3 or at the integer quantum-Hall transi- 
tion in d = 2, the states belong to the weakly multifractal 
regime, i.e. d — D2 <i d, the PBRM model allows for an 
investigation without such a limitation. The evolution of 
the generalized dimensions as a function of the parame- 
ter h therefore represent this behavior, i.e. Dq 1 for 
6^1 and in the other limit of 6 ^ 1 the multifractal 
dimensions vanish as Dq ~ bi^^ 

The multifractal dimensions, especially Di and D2 and 
their dependence on the parameter h in the case of the 
PBRM model, have been at the focus of several of our 
works li^iii In those works it has been demonstrated that 
a simple, phenomenological relation can be identified: 
Z?! w [1 + and D2 ~ [I + (a2&)"^]"^ where 

ai,2 are fitting constants. These continuous functions are 
trivial interpolations between the limiting cases of low-6 
and large-6, taking the half of the harmonic mean of the 
two cases. In Ref. [13 we generalized and propose the 
following heuristic expression for an extended range of 
the parameter q 



Dq^ [l + {aqh)-^]' 



(10) 



as a global fit for the multifractal dimensions Dq of the 
PBRM model for both symmetry classes, /3 = 1 and 
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FIG. 1: Dq as a function of 6 for the PBRM model at criti- 
cality with (a) /3 = 1 and (b) /3 = 2. The dashed Unes are fits 
of the numerical data with Eq. pup . 



(3 = 2. In Fig. [T] we show fits of Eq. ([T0| to numeri- 
cally obtained -Dg's as a function of b for some values of q 
and in Fig. [2ja) we plot the values of aq extracted from 
the fittings. We observe that Eq. ()10p fits reasonably 



well the numerical Dq for q > 1/2. It is important to 
stress that Eq. PH)) reproduces well the fe-dependencies 
predicted analyticall}*^ for the limits <C 1 and 1. 

The multifractal dimensions of Fig. [T] were extracted 
from the linear fit of the logarithm of the inverse mean 
eigenfunction participation numbers versus the logarithm 
of N, see Eq. ([T]). Di was extracted from the linear fit 
of the mean eigenfunction entropy versus the logarithm 
of TV, see Eq. We used system sizes of = 2" 

with 8 < n < 13. The average was performed over 2"""^ 
eigenvectors with eigenvalues around the band center of 
2i6-n ]-ealizations of the random matrices. 

By inspecting the obtained fittings using Eq. PH)) one 
may notice that Eq. ((6)) leads to 



X « (1 + aib)' 



(11) 



which also reproduces qualitatively well the b- 
dependencies predicted analyticallji^ii in the small- 
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FIG. 2: (Color online) (a) a, and (b) 7, = ai/a, as a function 
of q for the PBRM model at criticality with /3 = 1 and P = 2. 
The red dashed line in (b) equal to q is plotted to guide the 
eye. The error bars in (a) are the rms error of the fittings. 



and large-6 limits: 
1 -45 , 



X 



fi = l 
2 



6< 1 

6> 1 
(12) 

As a consequence of equating b in Eqs. ([T0|) and ([TT|) we 
get a direct relation 



2/37r& 



X « (1 - Dq) [1 + (7, - l)Dq 



(13) 



with 7q = ai/aq. We observed that jq ^ q in the range 
0.8 < q < 2.5, as plotted in Fig. [^b), so in this range of 
q values we can write simplified relations between x f^nd 
D„: 



l-Dr, 



l + {q- l)Dq 



and D„ 



1-X 



l + {q-l)x 



(14) 



The expression for Dq in Eq. (|T4)) reproduces Eq. ([7]) 
exactly for q = 1 and q ^ 2 and approximately for 1 < 
q < 2.5. Moreover, Eq. ([H]) combined with Eq. © allows 
us to express any Dq in terms of, for example, Di: 



Dq^Di [q+{l-q)DiY 



(15) 



We also noticed that by equating x for different Dg's 
from Eq. ()14p we could get recursive relations among 
them: 

q'Dq,{l-Dq,)-' =qDq{l^Dq)-\ (16) 

which in case of taking q' = q + 1 leads to 

Dq+, ^ qDqil + q - Dq)-' . (17) 
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Notice that all these relations can be expressed using 
the fact that the ratio qDq/{l — Dg) is independent of 
q. These expressions also provide a relation between the 
correlation dimension and the information dimension or 
between the correlation dimension and the compressibil- 
ity of the spectrum: 



(18) 



It is relevant to add that in the weak multifractal regime, 
i.e. when % — > 1, Eq. ([TS]) reproduces the relation given 
in Eq. (g]) with d = 1, reported in 



A. The case g < 1/2 

For q < 1/2, Eq. (fTO|) cannot be directly applied. How- 
ever, the regime q < 1/2 could also be explored within 
our approach by the combination of Eq. pU)) and the 
symmetry relationiS 

A, ==Ai_q with Ag = Dg{q-1) -d{q-l) , (19) 

implying that it is possible to link the multifractal di- 
mensions with indexes q < 1/2 to those with q > 1/2. 
Therefore, we get 



l-2q 



q 



ai-nh 



I — q \l + ai-qb 



(20) 



for q < 1/2; that is, once we know the coefficients aq 
for q > 1/2 we can use them to get Dq for q < 1/2. 
Moreover, by the use of Eq. (IT7|) it is possible to write 
Dq, for q < 1/2, as a function of any specific Dq with 
q > 1/2. For example. 



D„ 



Di 



l-2g q 
1-q l-q\l + q{Di-l) 



(21) 



provides Dq for q < 1/2 in terms of the information di- 
mension. Moreover, we can write down relations between 
X and Dq with q < 1/2: 



l~Dq 
q{2 - Dq) 



and 



Dr, 



1 - 9X 



(22) 



Finally, it is important to stress that Eqs. ()20II21|) get 
the form 



Dq = {l-2q){l-qr^ 



(23) 



in the limit 6 — > 0, which has been derived analytically 
in Refs. [Hi-ill . 



research activity for many decades Moreover, since 
the properties of the closed system, i.e. the fractality 
of the eigenstates, strongly influence the scattering and 
transport properties of the corresponding open system, 
the interest has also been extended to critical scattering 
systems. In particular, much attention has been focused 
on the probability distribution functions of the resonance 
widths and Wigncr delay timeS j^^'^^i^^" — as well as the 
transmission or dimensionless conductancei^"— 

Among many relevant results, here we want to focus 
on Wigner delay times t„ and recall that 

(i) for disordered systems at criticality the inverse mo- 
ments of Wigner delay times (tw'') scale a a^^'^^i^^ 

(t,;;«) cx N-"" where a, = qDq+i ; (24) 

(ii) for the PBRM model at criticality the typical values 
of the Wigncr delay times scale asii 



oc N" 



where CTi- = -Di 



(25) 



with 



= exp(lnTw). 

Note that Eqs. and (|25p provide a way to probe 
the properties of a critical system (i.e. the fractality of its 
eigenstates) by means of scattering experiments. More- 
over, we can also relate spectral properties to scattering 
properties by 



(i) combining Eqs. dH]) and (Pi)) : 
(T„ K, — where x 



(26) 



1 + 9X " 9(0-9 + 1) 

(ii) combining Eqs. (HH) and (^5]) : 

ar = l-x- (27) 

Also, we can express any Uq as a function of, say, ai: 

aq^2qai[{l + q) + {l-q)<j^]-K (28) 

Finally, note that we can obtain recursive relations for 
(Tg's in analogy to Eq. 



{q' + l)(jq,{q' -aq>) 1 = {q + l)aq{q~aq)- 



(29) 



or 



B. Wigner delay times 



<jq, = q'{q + l)aq[q{q' + 1) + (<? - q)aq]~^ , (30) 
which leads to 

aq+i = {q + lfaq[q{q + 2)-aq]-\ (31) 
when q' = q + 1. 



III. NUMERICAL RESULTS FOR THE PBRM 
MODEL 



As mentioned in the Introduction, the modeling and 
analysis of multifractal states in disordered systems at 
the Anderson-transition has been a subject of intensive 



In this section we present numerical justification of 
our analytical relations derived above using the PBRM 
model. 
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FIG. 3: (Color online) (1 - + (g - [see Eq. (HHl] 

as a function of b for the PBRM model at criticality with (a) 
/3 = 1 and (b) P — 2. The red dashed lines are the analytical 
predictions for x given by Eq. (|12|) . The blue symbols in 
(a) are independent numerically obtained values of x (taken 
from Ref. [48]'). Insets: qDq{l — Dq)~^ as a function of 6; see 
Eq. (|32|) . The red dashed line equal to oi\b is plotted to guide 
the eye. (a) ai = 4 and (b) ai — 9.4 were used. 




FIG. 4: Dq as a function of q for the PBRM model at criti- 
cality with (a) /3 = 1 and (b) j3 = 2. Black [Red] dashed lines 
are Eq. ^ [Eq. ([21])]. We used (a) [(b)] Di = 0.109, 0.429, 
0.843, 0.958, and 0.991 [Di = 0.230, 0.711, 0.931, 0.980, and 
0.996]. 



B. Wigner delay times 



A. Multifractal exponents 

In Fig.[3|we plot (1 - D,)[l + {q- l)Dq]-^ as a func- 
tion of h for several values of q for the PBRM model at 
criticality with /3 = 1 and {5 = 2 and observe good cor- 
respondence with the analytical prediction for x! that 
is, we verify the validity of Eq. (dH). For completeness 
in Fig. [SJa) we also include independent numerically ob- 
tained values of x (taken from Ref. [11]). In the insets 
of Fig. [21 we plot qDq{l — Dq)~^ as a function of 6, see 
Eq. pB|) , which for the PBRM model acquires the simple 
form 

qDq{l - Dq)-^ ^ aib . (32) 

The fact that all curves qDq (1 — Dq)~^ vs. b fall one on 
top of the other makes evident its independence of q. 

In Fig. [4|we present Dq as a function of q for the PBRM 
model with /3 = 1 and /3 = 2 for some values of b. As 
black and red dashed lines we also include Eqs. p3|) . for 
q > 1/2, and ([?T|) . for q < 1/2, respectively. In both 
equations we used values of Di we have obtained numer- 
ically. We observe very good correspondence between 
the numerical data and Eqs. ([^T|) and (jlSp mainly for 
-4 < q < 0.2 and 0.8 < q < 4, respectively. 



We obtain Wigner delay times Tw by turning the iso- 
lated system, represented by the PBRM model, into 
a scattering one by attaching one semi-infinite single- 
channel lead using perfect coupling. Since we are dealing 
here with the periodic version of the PBRM model, all 
sites are bulk sites and the place at which we attach the 
lead is irrelevant. To compute Tw we use the effective 
Hamiltonian approach described in Refs. [l7ll22| . For sta- 
tistical processing a large number of disorder realizations 
is used. Each disorder realization gives one value of Ty^ ■ 
We used N = 50, 100, 200, 400, and 800 getting 10^, 10^, 
10^, 10^, and 10^ values of r„, respectively. Then, the 
exponents aq [cr-r] were extracted from the linear fit of the 
logarithm of the averaged inverse moments of Wigner de- 
lay times [typical Wigner delay times] vs. the logarithm 
of N, see Eq. p4)) [Eq. ([25)) ]. We concentrate here on the 
PBRM model with 13 = 1 only. 

We start by noticing that if we combine Eqs. ([TU)) and 
()24p we get a heuristic expression for (jq as a function of 
b: 

'^1 ~ 1 I / ^ n-i • 

1 + (aq+ib) ^ 

In Fig.[S]we compare Eq. p3p to the numerically obtained 
CTq as a function of b for some values of q. We observe 
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FIG. 5: (Color online) (Jq/q (black open symbols) and -Dg+i 
(red full symbols) as a function of b for the PBRM model at 
criticality with P = 1- The dashed lines are Eq. (|33p : we used 
the values of a^+i reported in Fig. [2ja). For clarity the data 
for g < 3 was displaced upwards. 



that Eq. ([55]) fits reasonably well the numerical aq for q > 
0.1. In Fig. [5] we also include independent numerically 
obtained values of Dg that further verifies the validity of 
relation ^ [19,22.23]. 

In Fig. [5] we plot {q — (Tg)[q{aq + l)]^^ as a function of 6 
for the PBRM model at criticality; that is, we verify the 
validity of Eq. ([26| . We also plot the analytical prediction 
for X given in Eq. (jl2p and observe good correspondence 
with the numerical data. In the inset wc plot <Tq/ {q~ Uq) 
as a function of 6, see Eq. (f^ . which for the PBRM 
model acquires the simple form 



■"*-t'*^^ _ ' ' "o"'q = o'.l 
_i T, □ q = 0.5 

^1-%) \i \ O q=l 
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FIG. 6: (Color online) (g - cr,)[g((7, + 1)]-^ [see Eq. JMl)] 
as a function of 6 for the PBRM model at criticality with 
/3 = 1. The red dashed lines are the analytical prediction for 
X given by Eq. H12|) . Inset: aq{q — aq)~^ as a function of b; see 
Eq. (1311). The red dashed line equal to b is plotted to guide 
the eye. 
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(34) 



FIG. 7: (Color online) (Jr as a function of b for the PBRM 
model at criticality with (3 — 1. We include the numerically 
obtained Di, see Eq. (|25|) . The red dashed lines are 1 — X) see 
Eq. 1271). 



Finally, in Fig. |7] we show a-r as a function of b for 
the PBRM model at criticality with (3 = 1. To test the 
validity of Eq. (|27p we compare (Jt with the numerically 
obtained Di and with the theoretical prediction for 1 — x- 
We again observe good correspondence. 



IV. OTHER CRITICAL ENSEMBLES 

Remember that relations p4j|18p were obtained form 
the combination of Eqs. ([lUl) and pT|) . That is, rela- 
tions (|14m8p are expected to work in particular for the 
PBRM model at criticality. However, Eqs. repro- 
duce Eqs. ^ and which were shown to be valid for 
the PBRM model but also for other critical ensembles^ 
Therefore the question arises to which extent relations 
P4j|18p are valid for critical ensembles different to the 
PBRM model. So, in the following we verify the validity 
of Eqs. (|14m8|) for other critical ensembles. 



A. Calogero-Moser ensembles 

The Calogero-Moser (CM) A^-particlc systems yield 
three ensembles of N x N Hermitian matrices of the 
form^i^ 

Hmn = PmSmn + igi^ - 5mn)V{m - Tl) , (35) 

where Pm are independent Gaussian random variables 
with zero mean and unit variance, and V{m — n) is one 
of the three following functions 

1 1 1 

m — n ' sinh[(m — n)/A^] ' sin[(m — ri)/A^] 

These ensembles were denoted as^i^ CMR, CMH, and 
CMT, respectively. 

In Figs. |8| and |9| we plot Dq as a function of q for 
the CMR and CMT ensembles, respectively, for several 
values of g. To have an independent verification of our 
predictions, the data reported in these figures was taken 



7 




FIG. 8: (Color online) Dq as a function of q for the CMR 
ensemble for several values of g. The data was taken from 
Ref. d. Black [Red] dashed lines are Eq. ^ [Eq. ((2l|]. We 
used Di = 0.911, 0.771, 0.565, 0.226, and 0.0216. 
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FIG. 9: (Color online) Dq as a function of q for the CMT 
ensemble for several values of g. The data was taken from 
Ref. d. Black [Red] dashed lines are Eq. ^ [Eq. ^]. We 
used Di = 0.987, 0.949, 0.836, 0.416 and 0.0502. 



from Ref. 0|. We compare the numerical data with our 
equations for Dq with q < 1/2, Eq. ([2T|) . and q > 1/2, 
Eq. using as input, values of Di obtained by the 

interpolation of the curves Dq vs. q. We observe for both 
ensembles that our predictions reproduce reasonably well 
the numerical data. 



B. The Ruijsenaars-Schneider Ensemble and 
Intermediate quantum maps 

The Ruijsenaars-Schneider Ensemble (RSE) proposed 
in [43 is defined as matrices of the form 



Now, in Fig. [To] we present Dq as a function of q for the 
RSE for several values of g. The data was also taken from 
Ref. 0. As for the CM ensembles, here we observe that 
our predictions reproduce reasonably well the numerical 
data for both Dq with g < 1/2 and q > 1/2. The values of 
_Di we used as input in Eqs. (PT|) and ([T5|) were obtained 
by the interpolation of the curves Dq vs. q. 




FIG. 10: (Color online) Dq as a function of q for the RSE for 
several values of g. The data was taken from Ref. [^. Black 
[Red] dashed lines are Eq. ^ [Eq. ((2T|]. We used Di = 0.99, 
0.914, 0.75, 0.512, 0.193, and 0.022. 



200- 



100- 



Q 






o q = 1 


-r^ 20 


_ □ q = 1.5 




q = 2 


Q 


Aq = 2.5 




V q = 3 


10 


X q = 3.5 


Q 




cr 



(b) 



10 15 



20 



H„ 



cxp(i$,; 



1 — exp(27rig) 



iV[l - exp(27ri(m - n + g)/N)] 



(36) 

where 1 < m < n, $,„ arc independent random phases 
distributed between and 27r, and g is a free parameter 
independent on N. 



FIG. 11: (Color online) qDq{l — Dq) ^ as a function of g 
for (a) the RSE and (b) IQM for several values of q. Red 
lines are (a) Eqs. (|40I41|) and (b) Eq. ((42|. The muhifractal 
dimensions Dq reported in this figure were computed using 



the same matrix sizes and ensemble realizations as 
PERM model. 



for the 
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In fact, in Ref. [l8| we have already successfully tested 
the predictions given by Eq. ([T^ on the RSE. There, we 
even formulated heuristic expressions for Dq-. 



Dq^[l-{g~lf] [l + {q-l){g-lf]-' (37) 



and 



[e-{g-kf] [e + {q-l){g-kf] ' . (38) 



for < g < 1 and |g — fc| ^ 1 with k > 2, respectively. 
We got the expressions above by substituting x 
and X ~ (5 - ky/k"^ [or Di « 1 - (5 - 1)^ and Di w 
1 — [g ~ fc)^/fc^], obtained analytically in Ref. Q, into 
Eq. ([HI) [or Eq.JM]. 

Also, in Ref. [1^ we tested some of our predictions for 
a variant of the RSE, introduced in [50| , with the name of 
intermediate quantum maps (IQM) model; see also jsij . 
In this model the parameter g of the RSE is replaced by 
cN/g with cN = ±1 mod g, being g the parameter of the 
IQM model. For the IQM model we substituted x ~ 1 /g 
or Di w 1 — 1/g, analytical expressions reported in [50|, 
into Eq. or ([T5|) . respectively, to get 



taken from [521. We also include the prediction for Dg 
given by Eq. (where Dq has been replaced by Dq/2) 
using Di = 1.7405 ± 0.0002 We observe that the 

prediction of Eq. (jl5|) is a reasonably good approximation 



for Dq in the interval < q < 1.2. 

In Fig. lT^ b) we plot Dq and Dq'^ as a function of g, for 
the 3d Anderson model at criticality, together with the 
prediction for Dq given by Eq. (fTS)) (where Dq has been 
replaced by Dq/3) using Di = 1.93 ±0.01 [H. The data 
for Dq and U^^^ was taken from [s^. The multifractal 
dimensions Dq'^ were extracted from the scaling of the 
typical participation numbers TJ^^ = exp (Inio) with the 



system size N, from the relation Xq^'^ ^ iV 

We observe that the prediction of Eq. is reason- 
ably good for Dq'P with < q < 4. In contrast, Eq. 
do not reproduce the numerical Dq when q > 1. We 
have also substituted Di = 1.97 ± 0.002 (obtained from 
the interpolation of the Dq data from [sj] ) into Eq. ([15]) 
but the resulting Dq curve is very similar to that with 
Di = 1.93 ± 0.01, so we do not show it in Fig. [H] 



Dq^{g-l)ig + q-ir' 



(39) 



Here we just want to add that by the use of Eqs. ([57)1 
and §SI, for the RSE, and Eq. for the IQM model, 
we can demonstrate the independence of qDq{\ — Dq)^^ 
on q (already shown for the PERM model in Fig. [3]). 
In fact, by substituting the above-mentioned expressions 
into Eq. p6l) we get 



qDq{l-Dq)-' ^{g~l) 



\-2 



1 



and 



qDq{l-Dq 



fc2(g-fc)-2_l 



(40) 



(41) 



for < g < 1 and 
for the RSE; and 



|g — fc| ^ 1 with fc > 2, respectively. 



qDq{l 



1 



(42) 



for the IQM model. Then, In Fig. [IT] we plot qDq{l - 
Dq)~^ for the RSE and the IQM model for several values 
of q. We also include the equations given above in red 
dashed lines. We observe a rather good correspondence 
between numerical data and Eqs. (|40ll42p . 




C. Higher dimensional models 

The generalization of Eqs. p4j|18|) to higher dimen- 
sional systems {d > 1) can be done if Dq is replaced by 
Dq/d in Eqs. (|14j|18p . Below we explore the applicability 
of our results to the Quantum Hall transition in d = 2 
and the Anderson transition in d = 3. 

In Fig. [T^ a) we plot Dq as a function of q for the 
Quantum Hall transition (QHT). The data for Dq was 



FIG. 12: (Color online) (a) Dq as a function of q for the QHT. 
The red dashed line is the prediction for Dq given by Eq. p5|) 
using Di = 1.7405 ± 0.0002 [H. The numerical data for Dq 
was taken from [s^]. (b)-Dg and D^q^ as a function of q for the 
3d Anderson model at criticality. The red dashed line is the 
jrediction for Dq given by Eq. pS)) using Di = 1.93 ± 0.01 
531. The numerical data for Dq and D^'q^ was taken from 
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V. APPLICABILITY TO A DETERMINISTIC 
MODEL 

In the previous section we have verified that relations 
(I14j|18p are vahd for critical random matrix ensembles in 
Id and also, to some extent, to higher dimensional models 
at criticality. The common feature in the systems used 
above is the presence of multifractal eigenstates. How- 
ever, note that not only disordered models produce them. 
It is well known that deterministic models having self- 
similar potentials also possess multifractal eigenstates, 
see for example SSUsoj . Below, we test the applicability 
of our expressions relating multifractal exponents now 
to the multifractal eigenstates of a tight-binding model 
having a self-similar potential. 



^QoodooGockfOoedeeeee'eeeeef 
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FIG. 13: Dq, computed from Eq. (|43|l . as a function of q for 
the Fibonacci lattice for several values of g. The dashed line 
is the prediction for Dq given by Eq. H15|) with Di calculated 
from Eq. (|44|l . 



A. Off— diagonal one dimensional Fibonacci lattice 

According to Fujiwara et al. [s^ the multifractal spec- 
trum of a one-dimensional Fibonacci sequence can be 
represented by the inflation rule r„+i = r„T„_i, where 
Ti = A and T2 = AB, so T3 = ABA and so on. In this 
case the Schrodinger equation 

has a multifractal solution at the bandcenter, E = 0. 
Then, by deflning the parameter g = tAs/tAA the gen- 
eralized dimensions of the eigenstates take the form^ 

D,^{;i\na)-\q-l)-^\n[\\g^)/X{g^'^)\ , (43) 

where cr = (\/5 + l)/2 is the golden mean, and X{x) ~ 
{2x)^^[{x+ lY + \/{x + lY + 4x2]. Hence, the informa- 
tion dimension and the correlation dimension read as 

= (31na)-i [\Ti\{g^)^g'\ng^\'{g')l\{g^)\ (44) 



and 



= (31n(T)-i [21nA(g2)-lnA(/)] , (45) 

respectively. 

Then, in Fig.[T3]we plot Dq, computed from Eq. (|43p . 
as a function of q for the Fibonacci lattice for several val- 
ues of g. The dashed line is the prediction for Dq given 
by Eq. (fT5|) with Di calculated from Eq. (|44| . Again, 
as for the eigenstates of disordered models, here we ob- 
serve that Eq. ([T5|) reproduces rather well the multifrac- 
tal dimensions of the eigenstates of the Fibonacci lattice, 
mainly in the range of 1 < < 4. 



VI. CONCLUSIONS 

In this paper we propose heuristic relations on one 
hand between the generalized multifractal dimensions, 
Dq and Dqi, for a relatively wide range of the parame- 
ter q, and on the other hand between these dimensions 
and the level compressibility x- As a result we find a 
general framework embracing an earlier resullii^ and a 
recent onci^ Our proposed relations have been corrobo- 
rated by numerical simulations on various random matrix 
ensembles, as well as on a deterministic model having a 
self-similar potential, whose eigenstates have multifrac- 
tal properties. Of course the analytical relations and the 
numerical simulations set limitations on the validity for 
a certain range of the parameter q. Therefore our results 
are obviously approximate. Hence they call for further 
theoretical as well as numerical investigations. 

Moreover, since our relations between the generalized 
dimensions and the level compressibility allowed us to 
state a clear link between the spectral and scattering 
properties of disordered systems at the metal-insulator 
transition, it may be interesting to explore the con- 
sequences of our results on the quantities characteriz- 
ing the dynamical properties of critical random matrix 
ensembles; which have been the focus of very recent 
investigations!^"— 

We believe that our results may find applications in 
several recently studied models characterized by mul- 
tifractal eigenstates; such as deterministic self-similar 
potentials quantum spin chains^ Dirac fermions in 
the presence of random magnetic fields;^ and other crit- 
ical random matrix ensembles 1^1^ 



Acknowledgments 

The authors are greatly indebted to V. Kravtsov for 
useful discussions. This work was partially supported by 
VIEP-BUAP (Grant MEBJ-EXC13-I), PIFCA (Grant 
BUAP-CA-169), the Alexander von Humboldt Founda- 
tion, and the Hungarian Research Fund (OTKA) grants 
K73361 and K75529. 



10 



^ M. Janssen, Phys. Rep. 295, 1 (1998). Int. J. Mod. Phys. 

B 8, 943 (1994). 
^ F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 

(2008). 

^ A. Richardella, P. Roushan, S. Mack, B. Zhou, D. A. 
Huse, D. D. Awschalom, and A. Yazdani, Science 327, 
665 (2010); S. Faez, A. Strybulevych, J. H. Page, A. La- 
gendijk, and B. A. van Tiggelen, Phys. Rev. Lett. 103, 
155703 (2009); K. Hashimoto, C. Sohrmann, J. Wiebe, T. 
Inaoka, F. Meier, Y. Hirayama, R. A. Romer, R. Wiesen- 
danger, and M. Morgenstern, Phys. Rev. Lett. 101, 256802 
(2008). 

* E. Bogomolny and O. Giraud, Phys. Rev. Lett. 106, 
044101 (2011). 

^ E. Bogomolny and O. Ghaud, Phys. Rev. E 84, 036212 

(2011) . 

® E. Bogomohiy and O. Ghaud, Phys. Rev. E 85, 046208 

(2012) . 

I. Ruslikin, A. Ossipov, and Y. V. Fyodorov, J. Stat. Mech. 
L03001 (2011). 

^ E. Cuevas and V. E. Kravtsov, Phys. Rev. B 76, 235119 
(2007). 

^ F. Wegner, Z. Pliys. B 36, 209 (1980). 
1° J. T. Chalker, V. E. Kravtsov, and I. V. Lerner, JETP 

Lett. 64, 386 (1996); R. Klesse and M. Metzler, Phys. Rev. 

Lett. 79, 721 (1997). 
" A. D. MirHn, Y. V. Fyodorov, F.-M. Dittes, J. Quezada, 

and T. H. Sehgman, Phys. Rev. E 54, 3221 (1996). 
^2 A. D. Mirhn, Phys. Rep. 326, 259 (2000). 

V. E. Kravtsov and K. A. Muttalib, Phys. Rev. Lett. 79, 

1913 (1997); V. E. Kravtsov and A. M. Tsvelik, Phys. Rev. 

B 62, 9888 (2000). 
^'^ E. Cuevas, M. Ortuno, V. Gasparian, and A. Perez- 

Garrido, Phys. Rev. Lett. 88, 016401 (2001). 
1^ I. Varga and D. Braun, Phys. Rev. B 61, R11859 (2000); 

I. Varga, ibid. 66, 094201 (2002). 

J. A. Mendez-Bermudez, T. Kottos, and D. Cohen, Phys. 
Rev. E 73, 036204 (2006). 

J. A. Mendez-Bermudez and I. Varga, Phys. Rev. B 74, 
125114 (2006). 

J. A. Mendez-Bermudez, A. Alcazar-Lopez, and I. Varga, 
Europhys. Lett. 98, 37006 (2012). 

A. D. Mirlin, Y. V. Fyodorov, A. Mildenberger, and F. 
Evers, Phys. Rev. Lett. 97, 046803 (2006). 
^° C. Monthus and T. Garel, J. Stat. Mech. P09015 (2010). 
V. E. Kravtsov, private communication. 
J. A. Mendez-Bermudez and T. Kottos, Phys. Rev. B 72, 
064108 (2005). 

A. Ossipov and Y. V. Fyodorov, Phys. Rev. B 71, 125133 
(2005). 

A. Ossipov, T. Kottos, and T. Geisel, Europhys. Lett. 62, 
719 (2003). 

25 Y. V. Fyodorov, JETP Letters 78, 250 (2003). 

'^^ C. Texier and A. Comtet, Phys. Rev. Lett. 82, 4220 (1999); 

A. Ossipov, T. Kottos, and T. Geisel, Phys. Rev. B 61, 

R11411 (2000). 

T. Kottos and M. Weiss, Phys. Rev. Lett. 89 056401, 
(2002). 

2* F. Steinbach, A. Ossipov, T. Kottos, and T. Geisel, Phys. 

Rev. Lett. 85, 4426, (2000). 
2^ M. Weiss, J. A. Mendez-Bermudez, and T. Kottos, Phys. 

Rev. B 73, 045103, (2006). 



J. A. Mendez-Bermudez and I. Varga, Phys. Rev. B 74, 

125114 (2006). 
^1 B. Shapiro, Phys. Rev. Lett. 65, 1510 (1990). 

P. Markos, Europhys. Lett. 26, 431 (1994). 

P. Markos, Phys. Rev. Lett. 83, 588 (1999). 
^'^ K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 78, 4083 (1997). 
^5 K. Slevin, T. Ohtsuki, and T. Kawarabayashi, Phys. Rev. 

Lett. 84, 3915 (2000). 

K. Slevin, P. Markos, and T. Ohtsuki, Phys. Rev. Lett. 86, 
3594 (2001). 

X. Wang, Q. Li, and C. M. K. Soukoulis, Phys. Rev. B 58, 
3576 (1998). 

M. Riihlander and C. M. Soukoulis, Physica B 296, 32 

(2001) . 

M. Riihlander, P. Markos, and C. M. Soukoulis, Phys. Rev. 
B 64, 172202, 212202 (2001). 

I. Travenec and P. Markos, Phys. Rev. B 65, 113109 

(2002) . 

L. Schweitzer and P. Markos, Phys. Rev. Lett. 95, 256805 
(2005); J. Phys. A: Math. Gen. 39, 3221 (2006). 
M. Jansen, M. Metzler, and M. R. Zirnbauer, Phys. Rev. 
B 59, 15836 (1999). 

K. Senouci and N. Zekri, Phys. Rev. B 66, 212201 (2002). 

C. Monthus and T. Garel, Phys. Rev. B 79, 205120 (2009). 
"5 C. Monthus and T. Garel, J. Stat. Mech. P07033 (2009). 
■'^ J. A. Mendez-Bermudez, V. A. Gopar, and I. Varga, Phys. 

Rev. B 82, 125106 (2010). 

A. J. Martmez-Mendoza, J. A. Mendez-Bermudez, and I. 

Varga, AIP. Conf. Proc. 1319, 41 (2010). 

K. Slevin, private communication. 
■'^ E. Bogomolny, O. Giraud, and C. Schmit, Phys. Rev. Lett. 

103, 054103 (2009). 
^° J. Martin, O. Giraud, and B. Georgeot, Phys. Rev. E 77, 

035201 (R) (2008). 
5^ J. Martin, I. Garcia-Mata, O. Giraud, and B. Georgeot, 

Phys. Rev. E 82, 046206 (2010). 
52 F. Evers, A. Mildenberger, and A. D. Mirlin, Phys. Rev. 

Lett. 101, 116803 (2008). 
5'^ A. Rodriguez, L. J. Vasquez, K. Slevin, and R. A. Romer, 

Phys. Rev. Lett. 105, 046403 (2010). 
5* L. J. Vasquez, High precision multifractal analysis in the 

3D Anderson model of localization, PhD Thesis, The Uni- 
versity of Warwick, 2010. 
55 T. Fujiwara, M. Kohmoto, and T. Tokihiro, Phys. Rev. B 

40, 7413 (1989). 
5'' M. Woloszyn and B. J. Spisak, Eur. Phys. J. B 85, 10 

(2012). 

5^ V. E. Kravtsov, A. Ossipov, and O. M. Yevtushenko, J. 

Phys. A: Math. Theor. 44, 305003 (2011). 
5* V. E. Kravtsov, O. M. Yevtushenko, P. Snajberk, and E. 

Cuevas, Phys. Rev. E 86, 021136 (2012). 
5® I. Garcia-Mata, J. Martin, O. Giraud, and B. Georgeot, 

Phys. Rev. E 86, 056215 (2012). 

Y. Y. Atas and E. Bogomolny, Phys. Rev. E 86, 021104 
(2012). 

" C. C. Chamon, C. Murdy, and X. Wen, Phys. Rev. Lett. 

77, 4194 (1996); X. Chen, B. Hsu, T. L. Hughes, and E. 

Fradkin, Phys. Rev. B 86, 134201 (2012). 
^2 Y. V. Fyodorov, A. Ossipov, and A. Rodriguez, J. Stat. 

Mech. L12001 (2009). 



